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Abstract 

Consider the all-time maximum of a Brownian motion with negative drift. Assume that 
this process is sampled at certain points in time, where the time between two consecutive 
points is rendered by an Erlang distribution with mean 1/oj. The family of Erlang dis- 
tributions covers the range between deterministic and exponential distributions. We show 
that the average convergence rate as u> — > oo for all such Erlangian sampled Brownian 
motions is 0(w -1 / 2 ), and that the constant involved in O ranges from — £(l/2)/\/2~7r for 
deterministic sampling to l/v2 for exponential sampling. The basic ingredients of our anal- 
ysis are a finite-series expression for the expected maximum, an asymptotic expansion of 
y^j_i(l — exp(27rij/fc))~ s , s £ R, as k — > oo using Euler-Maclaurin summation, and Fourier 
sampling of functions analytic in an open set containing the closed unit disk. 
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1 Introduction 

Let {Bp{t) : t > 0} be a Brownian motion with negative drift whose position at time t is given 
by 

B p (t) = -pt + W(t), /3>0, (1) 

with Br(0) = and {W(i) : t > 0} a Wiener process (standard Brownian motion). Since 
/3 is assumed to be positive, the Brownian motion will eventually drift towards — oo, and the 
all-time maximum Mp = sup tgR + Bp(t) is well defined. In fact, it is known that Mg follows an 
exponential distribution with rate 2/3, so that P(Mg > x) = e~ 2 @ x (see e.g. [9l Lemma 5.5]), 
and hence the expected all-time maximum is simply given by MMg = 1/2/3. 

We consider sampled versions of the Brownian motion, meaning that we observe the process 
only at time points to = 0, t±,t2, ■ ■ ■■ A crucial assumption we make is that the times between 
consecutive sampling points T n = t n — i n _i, n£N, are independent and identically distributed 
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(i.i.d.). Let T denote a generic random variable with T = T\ (here = denotes 'equal in distri- 
bution') and ET = The constant u represents the expected number of observations per 
unit of time, henceforth referred to as the sampling frequency. 

It is readily seen that a sampled version of the Brownian motion constitutes a random walk 
{Sp{n) : n G N} with Sp{0) = and 

Sp{n)=X 1 + ...+X n with Xi,X 2 ,... i.i.d., X\=Bp(T). (2) 

The fact that Brownian motion evolves in continuous space and time leads to great simpli- 
fications in determining its properties. In contrast, the random walks that we obtain after 
sampling, moving only at certain points in time, are objects that are much harder to study. 
Although it seems plausible that, as ui — > oo, the behavior of the random walk should be similar 
to that of Brownian motion, there are many effects to take into account for finite oj. Let the 
maximum of the random walk be denoted by Mp(uS) = sup n=0 i Bp{t n ). The sampling error 
Ap(u) = Mp — Mp(uj) then depends on the drift /3, the sampling frequency u>, and of course 
on the distribution of T. This paper deals with the expected maximum of the random walks 
and, in particular, its deviation EA^(w) from the expected maximum 1/2/3 of the underlying 
Brownian motion. This relatively simple characteristic already turns out to have an intriguing 
description. 

We assume that the times between sampling points are drawn from an Erlang distribution, 

so that T = -Efc(A) with Ef.(X) an Erlang distributed random variable consisting of k independent 
exponential phases, each with mean 1/A, and 

fc-i 

P(£ fc (A) <x) = l-Y -e- Xx (\x) n , x>0. (3) 

n=0 

The random variable E^(X) has mean k/X and variance k/X 2 . One reason for working under 
the assumption of Erlangian sampling is that T constant and T exponentially distributed are 
opposite extremes with regard to randomness as well as in the family of Erlang distributions 
(viz. with X = kco and k — > oo and k = 1, respectively). Another reason is that Erlangian 
sampling leads to a random walk of which the distribution of the all-time maximum allows for 
an explicit solution. This gives rise to a series expression for ~EMp{u) = ¥,Mp(oj;k) in which 
the k terms involve the k roots of P(<r) = [cr(l + p — a)] k — p k = in \a\ < 1 with p G (0, 1) 
given by k(l — p) 2 ui = 2f3 2 p. In this paper this series expression is analyzed for the case that 
uj — > oo, and a major result is that 

EMp(w; k) = ±- ^g= + O^" 1 ) (4) 

where O(oj~ l ) holds uniformly in k > 1 as u — > oo, and tp^ — > ((1/2) as k — > oo. More than fifty 
years after its appearance, N.G. de Bruijn's Asymptotic Methods in Analysis [6], in particular 
Sec. 1.2 on the O-symbol, and Sees. 3.6-10 on the Euler-Maclaurin summation, still provides 
us guidance in doing the asymptotic analysis as required for establishing a result like @. For 
a more detailed overview of our results we refer to Subsection 12.21 Other recent works that 
exploit the beneficial properties of Erlangian sampling are [2] for compound Poisson processes 
and [12] for Levy processes. 

2 Preliminaries and overview 

In this section we present some preliminary results in Subsection 12. II and an overview of the main 
results in Subsection 12.21 Subsection 12.11 starts with results for the special cases of equidistant 



2 



and exponential sampling, and then we derive a general expression for the expected all-time 
maximum for Erlangian sampling. 



2.1 Special cases of Erlangian sampling 

In the case of equidistant sampling the time between two consecutive sampling points is always 
1/uj. From the definition of Brownian motion it then immediately follows that 

Bpffl+Ni-Pu,- 1 ,*,- 1 ), 

where N(a, b) denotes a normally distributed random variable with mean a and variance b. We 
should thus consider the maximum of a random walk with normally distributed increments, 
referred to in the literature as the Gaussian random walk. The maximum of this random walk 
was studied in OUT]. In particular, [Til Thm. 2] yields, for ui > /3/(2- v /7r), 

, N 1 C(l/2) P P 2 ^ C(-l/2-r) /-/? 2 V 
This implies immediately that 

EA^) = -^E + 0(^ 1 ) (6) 

with — £(l/2)/\/27r ~ 0.5826. Results similar to ([B]), in slightly different settings, have been 
presented in [H Thm. 2] and Thm. 1]. A crucial difference is that our result ([5]) is obtained 
from the exact expression for EMg(cj), while the results in [U [7] are derived from considering 
the Brownian motion in a finite time interval, and estimating its maximum by Euler-Maclaurin 
summation. 

In the case of exponential sampling, we assume that the times between consecutive sampling 
points are independent and exponentially distributed with mean 1/uj. In this case we can prove 
that (see Lemma 12. II with k = 1) 

Bp(T) = Ei(j 2 ) - £0.(71), (7) 

where 

7i = -p + + 2lj, 72 = P + VP 2 + 2w. (8) 

The random walk for which the increments are distributed as the difference of two exponentials 
has been thoroughly studied in the literature. The maximum of this random walk is known to 
be equal in distribution to the stationary waiting time in a so-called M/M/l queue with arrival 
rate 71 and service rate 72, for which (see e.g. [3l p. 108]) 

F(Mf,(u) >x) = 2V (72 - 7l)x , x > 0. (9) 
72 

This implies that 

EMp(oj) = K—!—, (10) 
72 72 - 7i 



from which it readily follows that 



EA (3 H = -^= + 0(^ 1 ) (11) 
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with l/y/2 ~ 0.7071. A similar result was obtained in Thm. 3] for a Brownian motion in a 
finite time interval sampled at uniformly distributed points. 

We next set T = Ej,(kuj) with mean uo~ l and variance (/cw 2 ) -1 . Notice that random sampling 
(k = 1) and equidistant sampling (k = oo) can be seen as special cases. We first make the 
following observation. 

Lemma 2.1. IfT = Ek{k(jj) then 

B /3{T) = E k (j 2 ) - -Efc(7i)> (12) 

where 

7i = -13 + \//3 2 + 2ku, 72 = (3 + yjl3 2 + 2ku. (13) 
Proof. For /3s + s 2 /2 < kuj, 

ku 



koj - (/3s + s 2 /2) 

7i 72 
(71 - s)(72 + s)' 



(14) 



Hence, 



E(e -.W.,^_2L_j ^ , (15) 

which completes the proof by Levy's continuity theorem for Laplace transforms [TH]. □ 

From Lemma 12.11 we conclude that in order to study an Erlangian sampled Brownian mo- 
tion, we need to study a random walk with increments defined as the difference of two Erlang 
distributed random variables. As it happens, this random walk has been studied before, and 
an explicit solution for the distribution of Mg(cj) is available. In order to explain this, we need 
to make a small excursion into the world of fluctuation theory. We start from the observation 

that Mp(uS) = W with 

W= max{0, W + B P (T)}. (16) 

This equation is a special case of Lindley's equation, describing the stationary waiting time of 
a customer in the GI/G/1 queue. The case in (|16p describes a single-server queue with Erlang 
distributed service times and Erlang distributed inter arrival times. In [1] it was shown that the 
distribution of W can be expressed as a finite sum of exponentials, the exponents of which are 
the roots of an equation that involves the Laplace transform of the Erlang distribution. This 
gives the result [H Corollary 3.3] 

fe-l 

F(Mp{u)>t) = J2 c 3 (J j e ~~ t2{1 ~ aj)t i t>0. (17) 

j=0 



Here, 



J YtiZli^i/o-j - 1) 

and o"o, . . . , Cfc-i the k roots in |<r| < 1 of 



, j = 0,1,..., k-1, (18) 



P(a):=[a(l + p-a)] k -p k = 0, (19) 
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where 

P = l-—, (20) 

72 

so that i = and 

fc(l-p) 2 w = 2^ 2 p- (21) 

The o"j are given explicitly as 



Oj = \ (1 + p) - ^/(| (1 + p)f - pe**V/k , j = 0,l,...,fe-l . (22) 
The expected all-time maximum then equals 

and it is this expression that forms the point of departure for this paper. Notice that (Tk-j = °j 
for j = 1, . . . , k — 1 and that <To = /) is real. Then it follows from (|18p that Cfc_j = for 
j = 1, . . . , k — 1 and that Co is real. Now (|23|) implies that EM^w; k) is real. 

2.2 Overview of the results 

In the coordinates p, k, u and /3, related according to (fT3|) . ([20]) . and (f2T]l . we obtain a limit 
result for EM^(w;fc) in the case that to — > oo and A; bounded or unbounded, and in the case 
that uj — > and k — > oo. In this paper, we do not address the intriguing question what happens 
when to tends to a non-zero finite limit and k — > oo. 

In Section [3] we start from (j23j) and show that when to — > oo and k is bounded, 

EM ^ k)= h-^( k -% 7t=) + 0( ^> <24) 



u . = e *rfj/* ) j = 0,1,..., A: - 1 (25) 



with 

the /c unit roots. Notice that the case k = 1 is in line with (jlip . 
In Section H] we determine the asymptotic behavior of the series 



fc-i 

yVl-u,)- (26) 



when s G R is fixed and fc — )■ oo. We use for this a method developed by Brauchart et al. [5] 
based on the Euler-Maclaurin summation formula. In particular, our results imply that 



k - y ■ ) = - * m «; 1/2) 1 + ( t -') . (2 7 ) 



fc-i 

\/2fc V ^ a/1 - Uj) V2tt 2V2 k 

Notice that for k — > 00 the right-hand side of (|27p tends to the leading-order term at the 
right-hand side of ([6]). 

In Section [3] the result (|24p is proved for the case that A remains bounded (Theorem 13. ip . 
The proof of Theorem 13.11 is relatively simple and uses the direct representation in (|18p of 
the Cj's as they occur in the series expression (|23p . A key observation in this proof is that 
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the term with j = in the series representation dominates all other terms. The result 
comprises the quantity k — X^=i(l — u j)~ 1 ^ 2 f° r which the asymptotic behavior as k — > oo is 
given (Theorem l4.1|) in terms of the Riemann zeta function by (|27|) . using an approach based on 
Euler-Maclaurin summation. It is therefore a relevant question to ask whether (|24p also holds 
with 0{oj~ l ) holding uniformly in k > 1. The approach to prove (|24p for unbounded k using 
the direct representation (fTBj) of the Cj is severely complicated by the fact that, as k — > oo and 
to — > oo, the k zeros of P in (|19p inside the unit circle accumulate on the set {l—y/l — u \ \u\ = 1}, 
so that the quantities o~i/o~j — 1, I ^ j, that occur in (|18j) . can become arbitrarily small. In 
a situation like this it may be advantageous, as exemplified on several occasions by N.G. de 
Bruijn in [6], to view the problem at hand from a different perspective. In this spirit, we use a 
different representation of the Cj's, viz. one in which factors crT" /a^ — 1 (rather than aj~ /crj — 1) 
appear with a~ and a + the zeros of P inside and outside the open unit disk \a\ < 1, respectively 
(Lemma I5.ip . This new representation gives bounds on c,- (Theorems 15.31 and 15. 4[) from which 
we can conclude that also in the case of unbounded k, the term with j = in the representation 
(|23p is dominant when lu — > oo. The proofs of Theorems 15.31 and 15.41 rely heavily on Lemma 
15.21 that gives bounds on the products H/iTo 1 ^ ~~ a T l a t) that appear in the new representation 
of the Cj . In the proof of Lemma 15.21 a crucial role is played by the result that expresses a 
series X^=oM u .j)j with h(z) = Y^=od n z n analytic in an open set containing the closed unit 
disk, in terms of the decimated coefficients d s k, s = 0, 1, . . . (Fourier sampling). In proving the 
extension of (|24|) to unbounded k (Theorem 16. ip , due to the results of Section attention can 
thus be restricted to the term with j = in (|23p . making the proof rather straightforward. 

In Section [7] we consider the behavior of EM«(w; k) as uj — > and k — > oo (the case that 
uj — > while k remains bounded yields EMg(u;;fc) — > in a trivial manner from (|18p - (|23p ). 
Then the term with j = in (|23|) is no longer dominant, and it can be shown from the results 
of Section [5] that EM^(w; k) is well approximated, see the proof of Theorem 17. 1\ by 

1 ~P V (1 + P ~ Qj) Qj 

2(3k ^(1-^)2(1 + ^-2^)- 1 j 

The series in (|28p can be cast into the form X^j=o -^(^j)' with F analytic in an open set 
containing the closed unit disk and P(0) = 0. To this series, the Fourier sampling technique, 
as it occurs in the proof of Lemma |5.2| can be applied. It thus follows that EM ( g(cj; k) tends to 
zero when oj — > and k — > oo, and k,uj related as in (|2ip with fixed (3 > 0, and also the rate at 
which EiW^u; k) tends to zero can be determined. 

3 Behavior of EMg(w; k) for bounded k and uj — > oo 

We prove in this section the following result. 
Theorem 3.1. Let k = 1,2, ... . As uj — > oo, 

k-l 



EM^;fc) = --^=(fc-^ Ji )+0(oj- 1 ), (29) 

)= 1 V 1 U 3 



1 1 

2/? ~ 



.r- 

where Uj are given in (|25j) . 

Proof. We use the series expression (f23|) for EM^(w;l:). We have from (fT3|) and ([20|l that 
1 — p = C^o;" 1 / 2 ), and so from 



1-a, = \(l- p ) + ^l-u ] -{l-p){l-u J ) + (\{l-p)Y 

= 71^7 + Oiuj" 1 ' 2 ) , j = 0,l,...,fc - 1 . (30) 
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In particular, 1 — do = O(oj 1 ^ 2 ) while 1 — aj is bounded away from for j = 1, k — 1 as 
ui — > oo. Then, from (fl~8]l and 1 — p = 0{uj~ 1 / 2 ), as oo — > oo 

C - \~ ^ 1 > ( 31 ) 



while for j = 1, k — 1 



IIKl^-i 



t^=- n n (^-0 < 32 > 

has a finite limit / as w -> 00. We conclude from p = do, ([3~T]) . ([32]) and 1 — p = 0(oj~ 1 I 2 ) 
that, as w — > 00, 

EM^u;; fc) = £ 7^ s = 1 ^ + O^- 1 ) • (33) 
72(1 - OjO 72 1 - P 

To proceed, we need to approximate Co accurately. Prom the identity in the first line of (|3U|) we 
find for j = 1, k — 1 

^ = l-i(l-p)- v /I^(l-i(l-p) + 0((l- y0 ) 2 )) 

= i-yr^77- + I( x /r^77--i)(i- p ) + O((i- /0 ) 2 ) . (34) 



fc-1 1 

1 - a-j 



This gives, see ((31 



pc = P k ] 



1 V T^--1( V T^--I)(l- P ) + 0((1- P r 



j TT V x "J 2 Vy ± J till ^^^Vv x ^ > 

= /(i + E4^)+o((i-p) 2 ) 

fe-i 

= l-(l-p)(k-Y,^^=)+0{{l-p) 2 ) , (35) 
3=1 V 1 n J 

where it has been used that p k = 1 — (1 — p)k + 0{{1 — p) 2 )- Using this in ([33]) while noting 
that 72(1 — p) = 2(5 and that ^- = , + ©(a;" 1 ), we obtain the result. □ 

The proof of Theorem 13.11 shows that the term with j = in the series (|23p dominates 
all other terms when u; — >• 00 and k remains bounded. We establish a similar result more 
generally, allowing k to be unbounded as well, in Section [61 and extend the result of Theorem 13. II 
accordingly. 

4 Asymptotics of X^=i (1 ~~ u i) s as k —>■ oo for fixed s 

The large-w expression in Theorem 13. II for EA/ J g(cj; fc) contains the series X^j=i (1 ~ u j)~ S with 
s = 1/2. It is of interest to find out how this series behaves with increasing k. Furthermore, 
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in Sectional we consider the case that u — > oo with unbounded k allowed, and then it appears 
that the behavior of the series Ylj=i (1 — u j)~ s f° r l ar g e k is required to be known for s = 
1/2,1,3/2,.... 

We adopt an approach in [5], for determining the asymptotic behavior of 

k-l 

S k (s):=^(l- Uj )- a , Uj = e 2 ^l k , (36) 
i=i 

as k — > oo and s £ M is fixed. In [5] this approach is used for finding the asymptotic behavior of 

k—l 

C4( S ):=]T(sin^)~ S (37) 

3=1 

as k — y oo and s £ C is fixed. The result that we obtain here for Sk(s) is of the same nature as 
the result for Uk(s) in [5], except that in our result, Theorem 14.11 below, all terms 

C(s-l)k s ~ l , 1 = 0,1,..., (38) 

occur, while the result in [5] has only terms (|38|) with even /. Furthermore, the exceptional 
cases s = 1, 2, ... are less complicated for our Sk(s) than they are for Uk(s) in [5]. 

The method in [5] is a fine application of the Euler-Maclaurin summation formula that can 
be found, along with various applications, in [SJ Sees. 3.6-10]. We take s & R in (f3"B"j) and this 
implies that the terms tj = (1 — Uj)~ s in ([36]) satisfy = t*, j = 1, ...,k — 1. Hence, Sfc(s) is 
real, and so we have 

fc-i 

S k (s) = £ / fe (j) , (39) 



where 
Expanding 



A0r) = Re[(l-e 2 ^'/ fc n . (40) 



— ) = £"1^' N< 27r ' ( 41 ) 
z=o 



with B\ the generalized Bernoulli polynomials, as in [5j Sec. 1], we get 



where 



f k (x) = J2 S > M< fc > ( 42 ) 

A(s)= (2ff) '~7' S)(0) + 1 = 0,1,.... (43) 



Remember that £ denotes the Riemann zeta function. 

Theorem 4.1. Let s £ K and Zei p = 0, 1, ... suc/i i/iai s + 2p > 0. Then, as k — >■ oo, 



2p 

fc + 2 ^ A(s) C(a " ^ + O^C**" 2 * -1 ) , « ^ 1, 2, ... , 

\k + 2 ]T AWC^-O^' + o.^"^ 1 ), « = 1,2,.. 



(44) 



where the constants implied by the O g)P depend on s and p but not on k. 



8 



Proof. We closely follow the approach in [Sj, so that many of the details are left out. We have 
for / G C 2p+1 ([l,m}) that 

E / + £(/(i) + /M) 

3=1 1 

p m 

+ E (tyr(/ {2 ^ 1) M-/ (2r " 1) (i)) + / f {2p+1 \x)^^d x , (45) 

r=l 1 

where C2 P +i(x) = B2 P +i(x — [x]) is the periodized Bernoulli polynomial and l?2r is the Bernoulli 
number. Using this with m = k — 1 and / = fk, see (|40|) . and noting that /fc(fc — sc) = fk(x), 
we get 

fc-i 

S*00 = E /fcW 

3=1 

fc/2 fc/2 

= 2 / AW^dHE^r 11 !^ / tf p+1) (-)§r|r da; - (46) 
i r - 1 i 

All quantities on the second line of (06]) involving fk can be expressed in terms of the fii{s) in 
In particular, we get 

fc/2 

/OO JLS— J 

f k {x) dx = k-2j2 #00 ^_ , (47) 



where it has been used that for s € K, s < 1 

1/2 





1 




^ Re 


(l _ e 2iviyy_ 





dy 1 (48) 



(1 _ e 27i%) 



and that the same analyticity considerations as in [5J Subsec. 2.1] apply. In (|47p we have to 
consider the cases that s = 1, 2, ... separately because of the term I = s — 1; this will be done 
below. 

We find, after using (|I2 >(P1) in (P) that for s / 1, 2, ... 

= * + 2 E aw {— m + a + E t^t (- - 02,-i 

Z=0 r=l V 

fc/2 

-(s-l) 2p+1 f - l - S - 2p - 1 ^^_d X }, (49) 



where (a) n is Pochhammer's symbol. The expressions in {. . .} at the right-hand side of (|49p are 
identified in [5j Subsec. 2.2] as incomplete zeta functions Cy,p(t) with y = k/2 and t = s — I, for 
which 

\Cy, P (t) ~ C(t)\ = OUy- 1 - 211 ) , y > . (50) 
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Hence, for s 7^ 1, 2, ... and any p = 0, 1, ... 

00 

5 fc ( S ) = fe + 2 £ C fc/2 > - k s ~ l . (51) 

In the case that s = K = 1, 2, ... , we take the limit s — )• K in (f5Tj) . using the result, to be proved 
below, 

„ fttiM j. (52) 

Note that Pk-i(K) = 0, due to the factor cos(^ + s)f at the right-hand side of (|43p . Hence, we 
get for s = K = 1, 2, ... 

00 

s fc (#) = \k + 2 ^ AW C k /2 >P ( K ~ • (53) 

(=0,(^-1 



Finally, taking any p = 0, 1, ... with s + 2p > 0, we can use (|50p to conclude the proof in the 
same way as the proof for Uk(s) in ([57|) is concluded in [SI Section 4]. 
We still have to show (|52p . From ()43p we have for if = 1, 2, ... 



s 



where5}f_ j 1 (0)/(K-l)! is the coefficient of z K ^ in (z/(e z - 1)) K . Thus, we have for < r < 2tt 

dz 



^S(O) 1 /■ 1 / * \* 



(K — 1)1 2m J z K \e z - 1 

\z\=r 

1 f / 1 \ K . 1 f 1 



2iri J V e z — 1 / 27rz y ur 1 + w 

\z\=r C, 



dz = — I -^—dw = (-l) K ~ l 



(55) 



where the substitution u) = e z — 1, = dw/(l + w) has been used and C r is the image under 
the mapping \z\ = r 1— > e z — 1 (which is easily seen, for small r > 0, to be a Jordan curve having 
the origin w = in its interior). This shows (|52p . □ 

Note 4.2. The result in (|27p can 6e obtained by taking s = 1/2, p = 1 in (|44|) . using that 

1 /tF 7T 3 / 2 

^ (l/2) = — , ^(1/2) = ^, A (i/ 2 ) = -_. (56) 

Note 4.3. /n i/ie lower case in (|44p a simplification occurs since for s = 1,2,... 

A(s) = 0, /c-Zodd; (( s -/) = 0, / = s + 2, s + 4, ... . (57) 

This leads to 

S k (s) = ^k + Yl &-ar( fl ) C(2r) A; 2r + O^Cfc-*" 1 ) . (58) 

r=0 

A similar situation as in Note 1*01 occurs in [5l Remark 1.2] for the case of Uk(s) in (|37p with 
s = 2,4, ... . It is concluded in [5] that one gets exact formulas for Uk(s) in that case. While 
this is true, see [10], the argument in [5] is incomplete. In the case of Sk(s) with s = 1, 2, . . ., it 
can be shown that it depends polynomially on k (degree < s), and so (|58|> holds exactly with 
the SiP deleted. 
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5 Bounds on cj from a representation using outer zeros 

When k is allowed to be unbounded, the analysis of EM«(w; k) using the series in ([23]) with the 
Cj given by (|18p is awkward. In this section, we present an alternative representation of the Cj, 
using the zeros of P in (|19|) outside |<r| < 1, that is crucial for the developments in this paper. 
For j = 0, 1, k — 1 , we let 



of = \ (1 + p) ± yj§ (1 + p)Y -pe^/k (59) 
be the two solutions of the equation 

a(l + p -a) =pe 2 ^l k . (60) 

Then for j = 1, k — 1 

\oJ\ <°o =P< 1 = °o < k/l > ( 61 ) 

and 

a" + at = 1 + p ; <7"a+ = p e 2 ^ . (62) 

When we let 



, j = 0, 1, k - 1 ; a j+k = af, j = 0, 1, A; - 1 , (63) 



then <7j, j = 0, 1, 2k — 1 are the 2A: zeros of P in (|19D , 
Lemma 5.1. For j = 0, 1, k — 1, 



Proof. We have from (JTHJ) that 

We re-express the product in the denominator at the right-hand side of ([65]). There holds 

k-l fc-1 

P(a) = (a(l + p - a)) k - p k = (-l) k (a - a;) (a - *+) . (66) 

1=0 1=0 

Hence, for j = 0, 1, k — 1 from ([63]) 

fc-1 

P'{a j ) = {-lf fa - *i) IJ " ■ ( 67 ) 

l=0,l^j 1=0 

On the other hand, 

P'(a 3 ) = (a k (l + p- a) k - p fc )'fa) = k p k ~ ^ . (68) 

Hence, for j = 0, 1, k — 1 , 



TT (n „\ (-1)*^) 



(-l)fcfe^ 1 + P~ 2 °i (69) 

^(l + p-a,) nfjo 1 fa - a+) 
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Next, by the first item in ([55]) . aj — = — o"/)> an d so 

n, . , i, 1 + p - 2(7,' 1 

Using (EOl) in (E5D, we get 



(70) 



1 (aj/p^jl + p-aj) 



k-l k-l 



(71) 



J k (l-a j )(l+ P -2a j ) l 

Finally, take out k factors from the last product at the right-hand side of (|7ip . and use, see 

{arf/pf = (aja+/p) k = 1 (72) 
to obtain the result. □ 

We now analyze the product Ylt^o ~~ a i/ <T ~j) °f which Yli=o (1 ~~ * s * ne special case 
with j = 0. 

Lemma 5.2. For j = 0, 1, k — 1, 

k-l 



1=0 



where 



Proof. We have 



Now 



Re( 5j (p))<0, \ 9j (p)\ <-±ln(l 



4p 



fe-i 

no . 

Z=0 J 



c - 



(1+p) 2 - 

-fcln<7+ + ^ M "/ - a i)j ■ 



k-l 



(73) 



(74) 



(75) 



1=0 



ln(cr+ - 07) = In 
= In 



o 2 - pe 27Ti i/ k + \/a 2 - P e 2 ™ l / k 



1 - aj + \fa 2 



where u\ as defined in ([25]) and 



a = i (1 + p) , 1 — a,j = \ a? — p u 



Hence 
where 
Note that 



ln ( a t ~ a = h j( u i) 
hj(u) := In 1 — aj + \J 'a? — pu 



T . = JL= Ap = 1 



a 2 (l + p) 2 \l + p 



1 - p\2 



\u\ < a 2 j p . 



€ (0,1) 



(76) 

(77) 

(78) 
(79) 

(80) 
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Furthermore, 1 — a,- + y/ 'a 2 — pu has positive real part when \u\ < a 2 / p = r 1 , and so hj(u) is 
analytic in an open set containing the closed unit disk. There is the power series representation 



(81) 



n=0 



in which d n (j) are the power series coefficients of In [1 — aj + \J a 2 — z] . 
From all this we get 



k-l 



k-l 



(82) 



1=0 

where it has been used that 

k-l 



£ 

1=0 



^irinl/k 



1=0 s=0 

{k , n = 0, k, 2k, 
, otherwise . 



(83) 



Noting that do(j) = hj(0) = lncr^ , we then see from (fT5j) that 



k-l 



II (l - 4) = ex P (* E <Wi)p fcs ) == expfe(p)) . 
z=o ct j 8=1 



(84) 



We shall derive an integral representation, see ([§0]) . for the d n (j), n = 1,2, ... , from which 
the bounds for gj readily follow. We have by Cauchy's theorem for n = 1,2, ... 



d n (j) 



1 f ln(l + aj — \/ a 2 — z) 



2iri 



r n+l 



dz 



(85) 



\z\=r 



when < r < a 2 . We deform the integration contour \z\ = r so as to enclose the branch cut of 
\/a 2 — z from z = a 2 to z = +oo. Now, for x > a 2 , we have 



•^/a 2 - (x ± id) = =Fz \J , 



x — a* 



(86) 



and so we get for n = 1, 2, . 



^n(j) = — : / hi(l — aj — i y/ x — a 2 ) — ln(l — a,- + i \/ x — a 2 ) 
2iri 1 



dx 

r n+l 



(87) 



By partial integration, noting that the quantity in [. . .] at the right-hand side of (|87p vanishes 
at x = a 2 , we get 



d n (j) 



+ 



27m J 2\J x — a 2 — aj — iy 'x — a 2 1 - aj + iy/x — a 2 



dx 
x™ 



-(l- aj ) f 1 



2irn 



dx 



y/x - a 2 (l-aj) 2 + x-a 2 x n 
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for n = 1, 2, ... . Finally, setting u — x — a > 0, we get 



oo 

1 — 0<f f 



dnij) = -—J (1 - a,) 2 + u 2 ^ • (89) 



o 



By the substitutions t = u/a and v = tb- , the result ([89]) can be brought into the forms 



dn(j) 



j 

00 j 1/2 



vrraa 2 ™ 7 1 + i 2 (1 + t 2 ) 1 
o 

oo 

-1 /■ 1 d« 



irna 2n J I + v 2 (I + v 2 /bj) 7 



n=l,2,... , (90) 



o 

where we have set bj = (a/(l — ctj)) 2 , so that 

b 3 ~ (i + P )2 e v 1 * 

is in the right-half plane. 

We now show that Red n (j) < 0. We have for t 2 = x > and b G C, Re 6 > that 

Re (lT^)= Re (F^) >0 - (92) 

Then Red n (j) < follows from O and Re bj 1 > 0. From ([SID it is then seen that Regj(p) < 0, 
and the first item in (|74p is proved. 

Next, we have from RefeJ 1 > that |1 + v 2 /bj\ > 1 for all v > 0. Hence, from (|9U|) . 

K(j)l< 2^' * = 1 > 2 '-- ( 93 ) 



From (|84|) it is then seen that 

l«WK*E5E3B = -H 1 -(3n. («) 



i 2ksa 2ks 2 V Va 



and noting that a = i(l + p), this gives the second item in (f74"j) . □ 

The following result follows immediately from Lemmas 15.11 and 15.21 
Theorem 5.3. For j = 0, 1, k — 1, 



1 

^ ^ k 



1 + p - a 



j 



(l-a j )(l + p-2a j ) 



(95) 



and 

^l ^^U^ +O^)) (96) 

where r = 4p/(l + p) 2 and i/ie constant implied by O in ([9"o]) is bounded by 1. 
Another inequality for c,-, j = 1, — 1 , is the following one. 
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Theorem 5.4. For j = 1, 2, [I 



--: m, 



I Cj | — | Cfc_ j | < 



1/2 



(97) 



Proof. We have from <7fc-j = a* , j = 1, k — 1 , that Cfc_j = c*, j = 1, k — 1 , and this gives 
|cj| = |c fc _j|, j = 1, ...,m. 

From Lemma 15,11 with j = we have 



co 



1 fc-i 



Therefore, as o"o = p, 



fe-i 



1-p/l 



Furthermore, from Lemma 15, 21 for j = 1, 2, k — 1 

fc-i 



1 TT n a2\ 1/2 1 - P 1/2 



n ' 

z=o 



< 1 , 



and so from Lemma 15. II 



\Cj\ < 



1 + p - a j 



(l-a j )(l + p-2a j ) 
We shall show that for j = 1, 2, m 

1/2 



1 ~ P 1/2 



from which the result follows at once. We have with t = 2irj/k 6 (0, 7r] 



l-^ = i(l-p) + A /(i(l + /) )) 2 - /0 e 



l + p-2<rj=2J(i(l + p))2-pe« 



l + j0 -o . = i(i + p) + ^/(i( 1 + p ))2_ pe « . 



Now 



Re 



> o 



and so 



|1 - aj\ > |(|(1 + p)) 2 - P e u \ 1/2 , |1 + P - 2*j\ > 2|(|(1 + p)) 2 - pe li | 1/2 
while, as < p < 1, 



l+p-<Tj| < i(l + p) + J(i(l + p)) 2 + p< 1 + V2 



(98) 



(99) 



(100) 



(101) 



|l-<Tj| > |1 + P-2<7j| > 2 (4) V2 ' + < 1 + , (102) 



(103) 

(104) 
(105) 



(i(l + p)) 2 -pe'* < J$(l + p))2 + p, (106) 



(107) 



(108) 
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which establishes the third inequality in (|1U2|) . For the first two inequalities in (|102|) . we compute 



|(|(1 + p)f - p^f = (|(1 + p)f - 2p(|(l + p)f cost + p 2 
= (1(1 - / ,))4 + 2 p(I(l + p)) 2 (l- cost) 

= (|(1 - p)f + p(l + pf sin 2 ^ > (|(1 - p)f + p(l + p) 2 (t/^ 2 , 



(109) 



where the inequality sin a; > 2x/tt, < x < ir/2 has been used. Now for < y < 1 and 
< p < 1 



( ' (1 / ))4 + P(l + P) 2 > " P)) 4 + P(l + P) 2 = ((|(1 + p)f + P? > ihf , 



It 



and so we get 



t \ 2 



|(I(1 + p)) 2 -pe*f>(±) 4 (£) 
This yields the first two inequalities in (|1U2D . 



< t < tt . 



(110) 

(111) 
□ 



6 Extension of Theorem 13.1 



In this section we show the following extension of Theorem 13.11 
Theorem 6.1. For fixed f3 > 0, 

fc-i 



: 4 - (^r) 1/2 { k ~ £ 7T=) + °^~ x ) 



2/3 \2kuJ 
where O holds uniformly in k = 1, 2, . 



uj — > oo , (112) 



Proof. We first show that we can restrict attention to the term j = in the series (j23|) for 
EM^w; k). We have by ([20]) 



k-l 



E 



x 72(1 - a,) 2/3 ^ 



1-p y4 Cjq-j 



< 



/5 ^'i-^ 



(113) 



where m = [^fc] and where it has been used that Ck-j = c*, cr^-j = a* , j = 1, 2, fc — 1. Now 
Theorem E31 and (fTU2"j) give 

1 - P | Cj(jj / 1 - P n . „i/2 /- \ Vk f2k\ 1 / 2 



.7=1 J .7=1 



< C V 2 C (3 )(1 _ p)2 , = (V2 + 2)cf C(l) ^ 



where (|21|) has been used in the last step. Hence, 



EMp(u J -k) = ^ + cl /2 l3 0(u- 1 ) 



(114) 



(115) 



where the constant implied by O is bounded by 2(\/2 + 2)£(|). 
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We now bound and approximate cq with (|124p . (|130|) as a result. We have from Lemma 15.11 
with j = 



co 



k-l 



0"-i 



Furthermore, there is the approximation, see the appendix, 

1 



l-<Tj = y/l-Uj (1-e 1 



+ 



(1 ~ P) ; 



with e = i(l - p) + §(1 - /o) 2 , and 

\^l^T j \ = \^/l-u k ^\>2(j/k) l l 2 , J = 1,2,..., 
We furthermore have from (1251) that 



m . 



k-l 



fc— 1 

[ sin — = 1 , 



fc 



A; 



(1-P) 5 



see [5], (1.11) for the last identity. From all this we get 

fc— l 

q ' = nO-^- 7 =)+-v xA ^- 

Next, from (|118p and e = 0(1 — p), we see that 

1 



el 



V 1 - u i 



O 



1 



uniformly in j = 1, 2, fc — 1 . Hence 
fc-i 



co =n i-h 1 



U 1 



(i - p) ; 



Now by (fTT8j) and m = [±fc] 



J=l ♦ 3=1 J 

< 1 fc^l - P )3 2( i = i fc(i _ p) 3 ^2 = 0(-^j- 2 



Hence 



with 



c = c ( 1 + O 



1 



w3/2jfel/2 



CQ 



exp( 2 gln(l- £ (l - 1 

i=i v 



(116) 



(117) 



(118) 



(119) 



(120) 



(121) 



(122) 



(123) 
(124) 



(125) 
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We develop 



fc-i 1 fe-i 

E^-^-^MWEf 1 -^- 



3=1 
fe-i 



fe 1 19 fe 1 i 



From 



3=1 v J 3 = 1 



(126) 



(127) 



(see (|27p for the first item in (|127p and use Theorem 14. II or proceed directly for the second item 
in (|127p ) we have 

1 N2 



fc-1 
3=1 



V 1 _ % 



(128) 



Finally, from (|118p we have 
fc-l 



V I 1 



Using (TrjSHBSD in (gZ5D , we get 



fc-i 



co 



3=1 



+ 0(e 2 k) + 0(e 3 A: 3 / 2 ) 



(129) 



3=1 



fc-1 

1 " (! " P) E f 1 - -71—) + + °( w_3/2 ) 

V U J 



3=1 



(130) 



where we have used that 2e = 1 — p + 0((1 - p) 2 ) and that k(l - p) 2 = O(co 1 ), see ([2Tp . 
Using ipn]l and (fT24"]l in (fTT5i we get 



EM„( W ;fc) = |+0(^; 



2/? 



fc-i 

3=1 v 



1 1-P 

2/3 2/3 



fe-i 



i + pE 1 



3=1 
fc-1 



+ 0(uT 1 ) 

+ o(^- 1 ; 



2/3 2/3 v ^ ,/r^7j; T v 



i V 1 -^- 



(131) 



where (|127[) has been used in the last step to replace the p in front of ^ • by 1 at the expense 
of an error 0((1 - p) 2 vTfc) = ©(a;" 1 ). Since by flZTj) 



1/2 



1 



1 - P = _P_ 

2/3 v / 2kJ V2k^ 



O 



we get the result. 



(132) 
□ 
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7 Behavior of EM^w; k) as uj — > and k — > oo 

In this section we show the following result. 

Theorem 7.1. Assume that j3 > is fixed and that u — > and fc — > oo. T/ien EM^(u; fc) — )■ 0. 
Proof. By Theorem 15.31 

1 + p - (Tj 



_ 1 

Cj ~ k (1 - o-j)(l + p - 2aj) 
where for all k = 1,2, ... and all j = 0, 1, k — 1 

4p 



(1 + e j>k (p)) 



We have by (J2T 



and so 



fc(i-p)" 



,2 



oo 



r k < exp 



-fc(l-/o) S 



(133) 

(134) 

(135) 
(136) 



Now by (EUD and (fT33l) 



fc-i 



EM^(c;; fe) 



1-/3 ^ 1 (1 + P ~ a j) a j 



E 



P£p,fc(p) 



fc-1 



E 



2pk(l - p) 2 2/3k ^ (1 - ^) 2 (1 + p - 2<7j) 



(137) 



The second term on the right-hand side of (|137|) tends to by (|134|) - (|136p . As to the third term 
on the right-hand side of (|137p . we estimate 



1-p 



fc-l 

E 



(I + p - aj) <jj e jtk (p) 



< 



2pk (l-aj) 2 (l + p-2aj) 

(l + p-aj)aj 



(l-p)r k 



k\ E 



2(3k(l-r k ) ^ i \{l-a J ) 2 (l + p-2a J ) 



j 

.k m 



- m-r k )^j_ 2 m 1/2 ^-^) 2 

J 2fc \2/c/ 



(138) 



Here we have used (|102p . with m = [^k], and \aj\ < 1. By (|135j) and (|136p we have that 
r fc (l — p) \fk — > 0, and so also the third term at the right-hand side of (|137p tends to 0. 
We finally consider the first term, 



Ri. :- 



1-p 



fe-l 

E 



(l + p-aj)aj 



2 P k f^ (l-<Tj) 2 {l + P~2aj 
on the right-hand side of (|137p . We show below that 



(139) 



< R k < 



1 + P 1 

2/3(1 - p) ^(k - 3/4) 1 - r k 



k - 



(140) 
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From (IT351) and (fTMjl it then follows that also R k -> 0. 

To show (jl40|) . we follow the approach that was used to prove Lemma 15.21 and we let for 

\u\ < t~ 1 

F(u) = , (14 l) 

V ; (l-a(u)) 2 (l + p-2a(u)) ' V ; 

where 
Using 
we have 



a(u) = a — \J a 2 — z ; a = |(1 + p) , z = pu . (142) 
(1 + p - <t(«) = pu , (143) 



pu 



(l-(7(«)) 2 (l + p-2a(«)) 2(1 - a + v 7 ^ 2 ^) 2 v 7 ^ 2 ^ 

oo 



n=0 

By contour integration as in (|85p ~ (|88p . we have that 

If z dz 



(J 



2vri J 2(1 - a + Ja^~z) 2 yfa^~z z n+1 
\z\=r 



OO 



1 /" (1 - a) 2 + a 2 - x 1 dx 



2vr 7 ((1 - a) 2 -a 2 + x) 2 Jx^a 2 x n 

a 2 
oo 

1 [ + p 2 )-x 1 dx 



2vr 7 (x - p) 2 Vx"^ 2 ^" 



(145) 



for n = 0,1,.... Now go = 0, see (|144p . and so the last integral vanishes for n = 0. The 
integrand in this integral changes sign once, from positive to negative at x = \{1 + p 2 ) > a 2 , 
and l/x n is positive and strictly decreasing in x > a 2 when n = 1,2, ... . It follows that g n > 0, 
n = 1,2, ... . Also, we have 



l(l + p 2 )-x ^ i(l+p 



2 ^ a 2 



(x - p) 2 " (a 2 - p) 2 (1 - p) 2 ' 
and so we conclude that for n = 1, 2, ... 



1 

< ff n < — 



< 2 \ 2 N2 — = 7- To , x > a 2 , (146) 



1 1 + p 1 



oo 

-I 


1 






\Jx — a 2 


X" 


1 


oo 

/ Vt- 


dt 


a 2 " 


1 


1 t n 



(147) 



From -lj F^(0) = <7 n p n and r = p/a 2 , we then get 
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We return to (fTMj) . As in ([55|) , we have 



Since F(0) = 0, we obtain from (|148|) that 



- 2/3 7T (i-p) 2 £W 

s-l 1 

oo 

1 l + p f (r/t) k dt 



2vr/3 1 - p J 1 - (r/t) fc V^ 7 ! ' 



(150) 



i 



Then using that 



we obtain 



(r/t) k T k 1 

oo 

1 l + p T k f dt , . 

0<r ><2^t4t^J ¥VT=i- (152) 

1 

For the remaining integral, we use the substitution i = e s , s > and the inequality e s l 2 — e~ s / 2 > 
s, s > 0, and we get 



oo oo 



I -jjL= = [ ;" (fc " 3/4)S ds < [ s- 1 / 2 e-^> ds . (153) 

J t k y/t^l J J e s/2 _ e -s/2 J 
v 



The last integral in (fT53l) equals (ir/{k - 3/4) ) l l 2 , and using this in (|152p we get (|140p . The 
proof is complete. □ 

Note 7.2. From the estimates of the three terms at the right-hand side of (|147j) . it is seen that 

EM^fw; fc) = ©(r^-lnr^') 1 / 2 ) . (154) 

A Approximating 1 — cr ? 

We present approximations of 1 — o~j, j = 1,2, ...,k— 1, that were needed at several places when 
(1 — p) Vk is small. With Uj as in ([25]) . we have 



1 - °j '- 2 



\{l-p) + ^{\(l + p)) 2 -p Uj 



i(i _ p) + yr^- ^ _ (i _ p) + 4 V ^- ) . (155) 



(i-p) 2 v/ 2 

2 v 1 ~ w -i- v 1 ~ '"j ^ ~ ^ _ w "•" ; 

We have for j = 1, 2, 



|1 — Uk-j\ = \1 — Uj\ = 2sin( -j-} > -j- , (156) 
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and so 

(1 " P? 



. < |(1 -pfk 1 ' 2 , i = l,2,...,fc-l . (157) 

V 1 - u i 

We develop the square root on the last line in (|155p under the condition that 

l-R<\, \{l-pfk^<\. (158) 



Then we get 



= i(i _ p ) + yr^-_ i(i - P )yr^-+ Ki-p) 2 
-l(i-p) 2 ^^d-p) 3 -^^ + - 

= v /r ^7( 1 ~K 1 ~ ,/T~H- ) (1 ~ p + ^ (1 ~ /))2) 

, (1-p) 3 
16^/1 - l'/a i - itj / 

= v^^7(i-K 1 --T^)( 1 -^ + i( 1 -^ 2 ) + °( i 7^=)) • ( 159 ) 





-Uj 


(1- 


p) 4 


1 - 


Uj 


1 





Uj 
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